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^— > Abstract. Light beams carrying orbital angular momentum, such as Laguerre- 

Gaussian beams, give rise to the violation of the standard dipolar selection rules during 
^ the interaction with matter yielding, in general, an exchange of angular momentum 

O" 1 larger than ft per absorbed photon. By means of ab initio 3D numerical simulations, 

we investigate in detail the interaction of a hydrogen atom with intense Gaussian 

i-H and Laguerre-Gaussian light pulses. We analyze the dependence of the angular 

momentum exchange with the polarization, the orbital angular momentum, and the 
carrier-envelope phase of light, as well as with the relative position between the atom 

\Q and the light vortex. In addition, a quantum-trajectory approach based on the de 

Broglie-Bohm formulation of quantum mechanics is used to gain physical insight into 

IT) the absorption of angular momentum by the hydrogen atom. 
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1. Introduction 

In the past few years a great interest has been focused on helical light beams that 
are able to transport orbital angular momentum (OAM) and spin in its propagation 
direction [HE] such as Laguerre-Gaussian (LG) beams. In fact, the transfer of OAM to 
matter has been already reported experimentally in Bose-Einstein condensates [3j and 
ensambles of cold atoms [4J , and has been investigated for molecular systems [5] . OAM 
has also attracted attention from the standpoint of applications ranging from optical 
tweezers [6j to novel quantum information protocols [7j. 

At variance with previous works, where the OAM was coupled to the center 
of mass of an atomic ensamble, here we focus on the interaction of an OAM light 
pulse with a single atom allowing for electronic transitions to both bound states and 
continuum. In this context, we have recently derived novel selection rules for light- 
matter interactions [8] where more than one h unit of angular momentum per photon 
are exchanged. Here, by means of 3D numerical simulations of the interaction between 
a hydrogen atom with an ultrashort and ultraintense light pulse carrying OAM, we 
investigate the angular momentum exchange in different photoionization scenarios and 
its manipulation through the light polarization (linear or circular), the carrier-envelope 
phase of the ultrashort pulse, and the relative position between the atom and the light 
vortex. 

Furthermore, we numerically calculate the de Broglie-Bohm quantum trajecto- 
ries j9j [lOj [11] of the temporal evolution of the electron wavefunction to obtain physical 
insight into the electron dynamics. This formulation has already been applied in the 
strong field regime as a method to reduce the computational time to obtain the dynamics 
of a multielectron system [T2l [13] or to study the interaction of the hydrogen atom with 
an intense laser pulse in ID or 2D, obtaining the high harmonic generation spectrum [14J , 
the above-threshold ionization energies of the electron [15] , as well as an insight on the 
role of the quantum potential in photoionization [16]. However, to our knowledge, this 
is the first time that 3D photoionization is addressed with de Broglie-Bohm trajecto- 
ries. Our motivation in using these trajectories focuses mainly on illustrating how the 
electron absorbs the angular momentum due to the light polarization (photon spin) and 
due to its transverse profile (OAM). 

The article is organized as follows. In section [2] we introduce the physical system 
consisting of a hydrogen atom interacting with an ultrashort light pulse described in the 
LG basis. We proceed investigating different photoionization scenarios, starting with 



the standard linearly and circularly polarized Gaussian pulses, section 3.1, We then 



compare these results with the ones with LG pulses, section [3~2t with linear or circular 
polarization, where the combination of spin and OAM is studied. We then investigate 
the role of the carrier-envelope phase and the case where the LG vortex is not centered 



with respect to the atom position in sections |3.3| and |3.4[ respectively. In section |3.5[ 
we examine longer light pulses, where trapping of the electron state by the light pulse 
is observed. Finally, we present the conclusions in section [4| 
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Figure 1. (colour online) Sketch of the physical system consisting of a light pulse 
interacting with a hydrogen atom at the origin, where the time pulse dependence and 
its associated initial parameters are introduced. As an example, we plot the transverse 
profile of the vector potential Q for a linearly polarized LG mode with topological 
charge £ — 1. 



2. Physical model 

The system under investigation is sketched in figure [Tj a hydrogen atom interacts with 
an ultrashort and ultraintense light pulse carrying OAM. Being the nucleus a massive 
object, we can assume that the dynamics are given mainly by the electron excitation. 
The non-relativistic Schrodinger equation describing the dynamics of the electron in 
interaction with the light field reads: 
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being ip(r,t) the electron quantum state at time £, V(r) the Coulomb potential due to 
the interaction with the nucleus, m the electron mass, q the electron charge, p the linear 
momentum operator, and A^^t) the vector potential of the light electromagnetic field. 



2.1. Ultrashort light pulse 

We consider an ultra-short laser pulse carrying OAM, whose spatial profile will be 
expressed as an LG mode [HE]. Considering the wave propagating along the z direction, 
these LG modes read 
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where p = \Jx 2 + y 2 , (j) = arctan(y/x), w(z) = w Jl + (z/z ) 2 is the beam width, 
Wq is the beam waist at z = 0, z = k wl/2 is the Rayleigh range, k is the carrier 
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wave number, R(z) — z{\ + (zq/z) 2 ) is the phase-front radius, &g(z) = — (2p + \£\ + 
1) arctan(z/z ) is the Gouy phase, and L]f'(£) are the generalized Laguerre polynomials 

p{q) i^/ J (p-m)\{\£\+m)\m\ q ' { ' 

The two independent indices £ = 0,±1,±2,... and p = 0, 1, 2, . . . correspond to the 
topological charge and the number of nonaxial radial nodes of the mode. Note that for 
£ = p = we recover the standard Gaussian beam. As shown in Refs. [H [2], in the 
paraxial regime LG modes carry a discrete OAM of £h units per photon along their 
propagation direction. 

The temporal envelope of our light pulse, see figure [TJ is parametrized by a sine 
squared function, with frequency uj e = n /N cyc r ) where N cyc and r are the cycle number 
and the period of the carrier wave, respectively. Hence, the vector potential can be cast 
in the form 

A(r, t) =A sin 2 (°^{z + a ) - ou e t) e 1 ^"^* LG e , p (P, <M , , 

x (Oh(z — ct + 7rc/uj e + a ) — Oh (z — ct + a )) + c.c, 
being 9h{x) the Heaviside or step function, a the Bohr radius, c the speed of light, x 
the 'global' carrier-envelope phase, uj = ck = 2tt/t the frequency of the carrier wave, 
and A the vector amplitude of the wave. In the Coulomb gauge, the vector potential 
has a small polarization component in the propagation direction [lj, which has been 
neglected here. Note that the wavefront is helicoidal. 

The polarization in the xy plane is included in A = A (ae x + /3 e y ), where a and /3 
are complex numbers satisfying |a| 2 + |/3| 2 = 1. For example, for a circular polarization, 
a = 1 and /3 = is. Note that here we make use of the notation of [17J, where a left 
(right) circularly polarized beam carries angular momentum sTi with s = +1 (s = — 1) 
along the propagation direction. A linear polarization forming an angle (j) with the 
x axis corresponds to a = cos (/) and /3 = sin cp that is an equal superposition of the 
two circular polarizations. We would like to emphasize the fact that the electric field 
amplitude and the LG 'local' carrier-envelope phase that interacts with each part of 
the electron wavefunction depend on the transverse position, at variance with the plane 
wave case. Note that since our description of the electromagnetic field is in terms of 
the vector potential A, it includes the effects of both electric and magnetic fields of the 
incident laser pulse. 

2.2. Selection rules 

Hamiltonian ([!]) of the system can be split as % = Ho + Hi + Hu, where Ho is the free 
hydrogen Hamiltonian (kinetic and Coulomb terms), Hi = —q(p • A + A • p)/2m and 
Hn = q 2 A 2 /2m. The selection rules associated to Hi are (see Appendix A): 

| AL\ < \£\ + 1 < L { + L f , AM = ±(£ + 5), AL + £ is odd. (5) 

It comes clear from ^ that for £ = we recover the usual selection rules for plane waves. 
Furthermore, for £ ^ we can excite transitions with an angular momentum exchange 



Transferring orbital and spin angular momenta of light to atoms 5 

of more than one unit of ft. Selection rules ([5]), in photon terms, can be thought as the 
absorption of a photon carrying a total angular momentum £ + s in the propagation 
direction, where s accounts for the spin angular momentum. 

In the case of plane waves, %n is uniform in the transverse plane, yielding just 
a ponderomotive force but, interestingly, for £ ^ it can produce transitions between 
spherical harmonics. In fact, proceeding in a similar manner as with %i ) one obtains [8]: 

| AL\ < 2\£\ <L { + L f , AM = ±2£, AL is even. (6) 

Since selection rules ^ appear from a term with A 2 we expect them to be associated 
with two-photon processes, and, therefore will be less significant than those associated 
with selection rules ([5]), since the intensity near the vortex, where the atom is located, 
is relatively small. 

In the following section, where the complete dynamics of the system are numerically 
solved, these selection rules will be essential to understand the interaction process. We 
would like to remark that these selection rules apply to scenarios beyond photoionization, 
such as transitions between bound states. 

3. Numerical simulations 

Initially, before the interaction with the light pulse, we assume the electron to be in the 
hydrogen's ground state, i.e., ^o(^) = ^{ r -> t = 0) = 1/ ' J7ia^e~ r ^ aQ . We will first consider 
the atom to be centered in the light propagation axis, i.e., at p — 0, corresponding to 
the maximum amplitude of the Gaussian profile or the LG vortex position, and a pulse 
carrier-envelope phase X = 0. Eventually, we will assume arrangements where the atom 
is not centered in the light vortex and study the role of the carrier-envelope phase. 

In all the following cases we consider the pulse, see pi), to have N cyc = 3 (except 
in section 3.5), a carrier frequency of u = 1 au (2n x 6.57 • 10 15 s _1 , ultraviolet) and 
period r = 2n au (152 as)[§J We assume a beam waist that satisfies the paraxial regime, 
Wq = 9 • 10 4 au (4.79 /xm), which is much larger than the typical size of the atom 
(wq 3> a ). Note that in this regime, the Gaussian case is equivalent to a plane wave. 
For simplicity, since the OAM of the light is given only by the topological charge £ ) we 
will restrict ourselves to LG modes with p = 0. 

We express the electron quantum state at each time as \ift(t)) = c |^o) + \Si/j(t)), 
where \5i/j) is the excited part of the state and c = (ipo\ip(t)). Since the field has a 
frequency uj = 1 au, larger than the bound hydrogen energy (0.5 au), we assume that, 
at the end of the pulse, the excited part will mainly account for ionization. Thus, we 
define the ionization probability as Pi = \(5ip(t)\5ip(t))\ 2 = 1 — |c | 2 . To clearly observe 
the effects of the light on the atom, we will only plot in the figures the probability 
density of the excited part of the wavefunction. 

§ Atomic units (au), where ft = m = q = l/(47reo) = 1, are used throughout the paper. The length 
unit is ao and the unit of angular momentum is ft. 
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By means of the Crank-Nicolson algorithm, we integrate Schrodinger equation (fTl) 
in a 3D grid of N x x N y x N z = 400 x 400 x 250 points, covering a space of (—20, 20) au 
x (—20, 20) au in the polarization plane and of (—12.5, 12.5) au in the propagation 
direction. We take absorbing boundary conditions in order to avoid unphysical 
reflections when the ionized wavefunction reaches the borders of the integration grid. 
The decrease in the norm of the wavefunction due to these absorbing boundaries will 
be in most cases relatively small (less than 2 %) and as usually, it will be regarded as 
ionized population with kinetic energy large enough to reach the boundaries. 

During the evolution we will calculate average values for different observables such 
as the position of the electron (r), or its angular momentum (L) = (r x (p — qA)). 
Moreover, we will project numerically the quantum state onto the spherical harmonics 
in order to check the agreement with the previously presented selection rules. To acquire 
physical insight on the dynamics, we will compute the de Broglie-Bohm quantum 
trajectories from the evolution of the quantum state (see Appendix B for details of 
the quantum trajectories formulation). To ease the reading, from now on we will refer 
to them simply as trajectories. 

3. 1 . Gaussian pulse 

In this section, we consider a Gaussian pulse, i.e., with £ = 0, p = 0. We take the 
electric field amplitude to be A uj = 0.33 au (corresponding to a peak intensity of 
3.9 • 10 15 W/cm 2 ) for both cases of linear and circular polarization. 

3.1.1. Linear polarization The results of the simulation with a pulse polarized in the 
x direction, plotted in figure [2j show that the electron begins to be ionized in the first 
cycle and, at the end of the pulse, 30% of the quantum electron state is ionized, see 
figure [21(b). We also see in figure [21(c) that there is no exchange of angular momentum 
in the z direction since both the light pulse is linearly polarized and does not carry 
OAM. Figure \2Xd) shows that the most excited spherical harmonics are Y± and Yi _1 , 



in complete agreement with the selection rules discussed in section |2.2| since the linear 
polarization is a superposition of left and right circular polarizations. 

The motion of the electron wavepacket is unravelled in the trajectories plotted in 
figure[2](e). Trajectories describe oscillations at the field frequency along the polarization 
axis following (with a n/2 phase delay) the oscillation of the electric field. Some of the 
trajectories, i.e., the ones that start far enough from the nucleus, are ionized near the 
electric field maxima. The mean value of the position oscillates in the x axis while 
remains in the origin of the yz plane, except for a small shift (<~ 10 -2 au) on the z 
direction due to the magnetic field, in agreement with previous results [IB]. A shift 
around this value will appear in all the following scenarios. 

The Coulomb potential acting over the trajectories that are closer to the nucleus 
retains them in the nucleus vicinity. The ones that distance from it make the trajectory 
cloud expand. This can also be seen as the spreading of the quantum state function 
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Figure 2. (colour online) Results of a simulation with a Gaussian pulse (£ = 0, p = 0) 
linearly polarized in the x direction with electric field amplitude Aqu = 0.33 au. For 
the rest of the parameters, see text, (a) Projection of the excited state \Si(j(t)) onto the 
plane xy at times t — r, t = 2r, and t = St. Temporal evolution of (b) the ionization 
probability Pi and (c) the expected value of the angular momentum along the z axis, 
(L z ). (d) Projection of the electron quantum state at the end of the pulse into spherical 
harmonics Yjf . Arrows correspond to transitions allowed by selection rules shown in 
section |2.2| (e) Projection of the 3D de Broglie-Bohm trajectories onto the x and y 
axis over time and in the xy plane. The thick black curves in (e) correspond to the 
mean value of the electron position. 



away from the nucleus where the Coulomb potential vanishes, see figure [2](a). 



3.1.2. Circular polarization In figure [3] we present the results of a simulation with 
a right circularly polarized pulse. We have checked that in the case of left circular 
polarization the results of the simulation are completely symmetric to the ones discussed 
below. 

As in the case of linear polarization, the ionization probability reaches about 30% 
at the end of the pulse, see figure [3j At variance with the linearly polarized Gaussian 
pulse, which does not transfer angular momentum in the z direction to the electron, 
in figure |3](c) we see that, after the interaction, the electron has (L z ) = —0.3 au. As 
selection rules predict, we expect that a right (left) circularly polarized light mostly 
excites the Y^ 1 (Y±) spherical harmonic, which is in accordance with the numerical 
results in figure [3](d). Note that the electron is excited to L = 1 states through one- 
photon processes and to L = 2 states through two-photon processes. 

The trajectories depicted in figure [3](e) describe oscillations along the two 
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Figure 3. (colour online) The same as for figure^] but with right circular polarization 
(s = —1). The arrow in (e-iii) indicates the rotation direction of the position mean 
value. 



polarization axis. It is noteworthy that for each of the trajectories, the main dynamics in 
these two axis are independent. The mean value of the electron position shows that the 
electric field displaces the entire trajectory cloud from the propagation axis and makes it 
rotate clockwise (since it is right circularly polarized) around its initial position. This is 
due to the fact that the electric field is almost uniform in the transverse plane, therefore 
acting in the same way over all the trajectories. 

3.2. Laguerre-Gaussian pulse 

From now on we consider pulses bearing OAM, e.g., with £ = 1 and p = 0, for both the 
linear and circular polarization cases. 

We consider our atom to be centered at the light vortex, where the electric field 
amplitude is zero and increases linearly in the vicinity of the singularity. Thus, we 
require very intense lasers to affect the atom. We choose a pulse with A uj = 1.4- 10 4 au, 
corresponding to a peak intensity of 6.7 • 10 24 W/cm 2 , although the electric field at a 
distance of 1 au from the vortex is only 0.16 au. Albeit the considered light intensity 
is unreachable nowadays, the reader must take into account that very strong lasers are 
under construction [T9] . planning to reach up to 10 26 W/cm 2 at 800 nm wavelength. 



3.2.1. Linear polarization In figure [4] we show the results of a simulation with an LG 
pulse with £ = 1 polarized in the x direction. We have checked that due to the symmetry 



Transferring orbital and spin angular momenta of light to atoms 







Figure 4. (colour online) The same as for figure El but with an LG mode with £ = 1 
linearly polarized in the x direction with electric field amplitude Aqlj = 1.4 • 10 4 au. 



of the problem, the interaction with an £ = —1 pulse leads to the expected symmetric 
results. 

As shown in figure |3|(a), during the interaction with the pulse the electron is ionized, 
but it remains trapped around the light vortex, due to the ponderomotive potential of 
the pulse. Once the electric field intensity decreases, the electron state begins to expand, 
as it can be seen in figure |3|(a-iii) or the trajectories depicted in figure[I](e). This trapping 
effect will be studied in more detail in section |3.5[ Note that even though the pulse is 
linearly polarized in the x direction, the OAM breaks the symmetry on the y axis. At 
the end of the pulse the ionization probability is 53 %. 

It is interesting to point out the total angular momentum transferred to the electron, 
see figure [41(c). The electron starts in the ground state, with zero angular momentum 
and, as the pulse interacts with the electron, it oscillates, reaching, at the end of the 
pulse, a finite amount of angular momentum of 1.53 au, at variance with the linearly 
polarized Gaussian case where no exchange of angular momentum took place. Thus, we 
expect that the populated electron excited states at the end of the pulse bear angular 
momentum. The spherical harmonics populations plotted in figure [41(d) show perfect 



agreement with the selection rules from section |2.2| since the linear polarization is a 
superposition of the two circular polarizations (see next section for further details), i.e., 
AM — £ ± 1. Note that the electron is excited to L = 2 states through one-photon 
processes and to L = 4 states through two-photon processes. 

The mean motion of the electron state is negligible. Nevertheless, the dynamics 
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Figure 5. (colour online) The same as for figure Hmut with right circular polarization 
(* = -!)■ 



are far from simple, due to the electric field having a spatial profile where it points in 
opposite directions in opposite points of the polarization plane. A peculiar effect is that, 
even though the polarization is linear in the x direction, the trajectories are forced to 
acquire a velocity in the y direction, in order to rotate counterclockwise, due to the light 
£ = 1 OAM, giving rise to the absorption of angular momentum by the electron. 

Notice that in the circularly polarized Gaussian case, the rotation of the trajectories 
around the propagation axis took place in sequence with each other, whereas here each 
of the trajectories rotates independently around the axis. 



3.2.2. Circular polarization In this section we consider the same pulse as in 



section 3.2.1 (£ = 1), but circularly polarized. At variance with section 3.1.2, the 
dynamics for the left and right circular polarization cases are different due to the relative 
sign between the light OAM and polarization components of the angular momentum 
along the z direction. This comes clear from the results shown in figures [5](a) and[6](a), 
whereas the right circular polarization preserves the circular symmetry of the excited 
quantum state, the left circular polarization gives rise to a more complex behaviour 
due to the exchange of angular momentum. Note that in the case of right circular 
polarization (s = — 1) the spin and OAM angular components compensate while in the 
left circular polarization case (s = 1) they add up. 

Figures [5](b) and [6](b) show that, even though the electric field intensity is the 
same for both cases, in the left circular polarization case the ionization is approximately 
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Figure 6. (colour online) The same as for figure Hmut with left circular polarization 
(* = +!)■ 



twice that in the case of right circular polarization (63 % and 31 %, respectively). In 
this particular case, we attribute this enhancement of the ionization to the exchange 
of OAM, as can be seen comparing figures [5](c) and[6](c). At the end of the pulse the 
expected value for L z in the left circular polarization case is 3.35 au, while becomes null 
for the right circular polarization. 

The transfer of angular momentum to the electron is clearly seen in figures [51(d) 
and[6](d), in perfect agreement with the previously discussed selection rules. We can 
check that a right circular polarization excites states with M = 0, as AM = £ + s = 0, 
therefore not transfering any OAM in the z direction to the atom. On the other hand, 
a left circular polarization excites states obeying AM = AL = £ + s = 2. 

As in the linear polarization case, the main dynamics are concentrated on the 
polarization plane and the mean motion of the electron state is again negligible. The 
absorption of angular momentum by the electron is clearly depicted in the trajectories: 
while in the left circular polarization the electric field of the light exerts a torque on the 
trajectory cloud forcing it to rotate around the vortex, see figure [6](e-iii), the combination 
of OAM and right circular polarization partially inhibits this rotation see figure [5](e-iii). 



3.3. Carrier- envelope phase manipulation 



Here we consider the same pulse that in section 3.2. 1| (a linearly polarized £ = 1 LG 



pulse), but modifying the carrier-envelope phase (CEP) x, see expression of the vector 
potential Q. In figure [7] we show the results for two different values of the CEP (a) 
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Figure 7. (colour online) Results of a simulation with an LG pulse (£ = 1, p = 0) 
linearly polarized in the x direction with electric field amplitude Aquj = 1.4 • 10 4 au, 
for a carrier-envelope phase of (a) x — n an d (b) x = f ■ Temporal evolution of (i) the 
ionization probability Pi and (ii) the expected value of the angular momentum along 
the z axis, (L z ). (iii) Projection of the electron quantum state at the end of the pulse 
into spherical harmonics Yjf . 



X = 7T and (b) x = 7r /2- Now, the ionization reaches 52% for x = ^ an d 55% for 
X = 7r/2, very similar to the pulse of section 3.2.1, of 53 %. The fact that the ionization 



probability does not significantly depend on the CEP for linearly polarized LG pulses, 
can be understood by noticing that the 'local' CEP that interacts with each part of the 
electron wavefunction takes all possible values due to its azimuthal phase dependence. 
The final angular momentum is 1.30 au when x = ^ an d 1.39 au when x = 7r /2, 
relatively smaller than the 1.53 au obtained obtained for x — 0- The fact that the 
absorbed angular momentum by the electron slightly varies with the CEP affects the 
probability of each excited spherical harmonic, but not the states permitted by the 
selection rules, which are the same in all cases, as can be seen comparing figures [7](a-iii), 
[Tfb-iii) and UK 



3.4- Off-center vortex 



In this section we consider again a linearly polarized £ = 1 LG pulse, as in section [3.2.1 , 
but taking into account a small displacement of the light vortex, while keeping the 
hydrogen position at the origin. The results of a simulation where the vortex is located 
at (x, y) = (0, 1) au are plotted in figure[8](a). From the projected density of the electron 
wavefunction is difficult to perceive the transverse displacement of the electron, so we 
plot the mean value of x and y in figure [8](a-i) to study the motion of the electron. 
The expected value of the electron position shows that the motion begins oscillating in 
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Figure 8. (colour online) Results of a simulation with an LG pulse (£ = 1, p = 0) 
linearly polarized in the x direction with electric field amplitude Aquo = 1.4- 10 4 au, for 
the light vortex placed at (a) (x, y) = (0, 1) au and (b) (x,y) = (1,0) au. (i) Expected 
value of the electron in the xy plane, ((x), (y)). The cross indicates the position of the 
vortex. Temporal evolution of (ii) the ionization probability Pi and (iii) the expected 
value of the angular momentum along the z axis, (L z ). (iv) Projection of the electron 
quantum state at the end of the pulse into spherical harmonics Y^f . 



the polarization axis until it is far enough from the nucleus and then starts rotating 
counterclockwise around the vortex. In this case, at the end of the pulse, Pi = 62 % 
and (L z ) = 1.86 au, see figures |8|(a-ii) and|8](a-iii), both greater than in the case where 
the vortex is located at the origin, because now the atom starts at a position where the 
electric field is more intense. 

On the other hand, figure [8](b) shows the results for the light vortex located at 
(x,y) = (1,0) au. We can see a clear mean motion of the electron too, mostly in the 
polarization direction but the OAM breaks again the symmetry making the electron 
rotate counterclockwise around the vortex. As in the previous case, at the end of the 
pulse we obtain larger values for Pj (65 %) and (L z ) (1.79 au) than in the case where 
the vortex is located at the origin. The analysis of the spherical harmonic spectrum 
of the excited electron state is shown in figures J8](a-iv) and |8](b-iv) . Comparing these 
projections with the ones in the case of the optical vortex placed at the origin, see 
figure |4](d), we see the appearance of a small contribution of YJ 1 , Y^ 1 and Y 3 3 , showing 
the excitation of electronic transitions with AL = 1 that, at first glance, seem to be in 



disagreement with the selection rules. The selection rules from section [272] are calculated 
for the atom located at the center of the optical vortex. Thus, to address the problem 
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Figure 9. (colour online) Results of a simulation with an LG pulse (JL — 1, 
p = 0) 10 cycles long linearly polarized in the x direction with electric field amplitude 
Aou = 1.4 • 10 4 au. (a) Temporal envelope of the pulse, (b) Projection of the excited 
state \5ip(t)) onto the plane xy at times (i) t = 3r, (ii) t = St, (iii) t = 9r and (iv) 
£ = lOr (c) Temporal evolution of the standard deviations in the x (solid red curve) 
and y (dashed green curve) directions of the electron distribution |^(r, t)\ 2 . 



appropriately, we would need to rewrite the vector potential in terms of LG modes 
centered at the center of the atom in order to take into account the displacement of the 
light vortex. This would lead to components with £ = 0, and thus, to the appearance of 
the electronic transitions with AL = 1. Note that an analogous interpretation could be 
made by writing the electron state spherical harmonics centered at the vortex position. 
The study of the full dependence of the absorption of OAM on the relative position 
between atom and light vortex is computationally very demanding and it is out of the 
scope of this paper. However, the small displacement of 1 au here considered has been 
enough to show that electronic excitations with large angular momentum take place 
even if the electron is not on the optical vortex. In any case, increasing the distance 
between the atom and the light vortex up to a point that the azimuthal dependence of 
the field becomes negligible, it is reasonable to expect that the standard selection rules 
for plane waves would be recovered. 



3.5. Longer pulses 



In this section, we consider the same linearly polarized £ = 1 LG pulse as in section 3.2.1 



but with a pulse 10 cycles long. Thus, the sine squared envelope in Q is substituted 



u env (t) 
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by: 

'sin 2 (^/(6r)) for t < 3r 

1 for 3r < t < 7t , (7) 

sin 2 (tt(10t - t)/(6r)) for £ > 7r 

i.e., the temporal envelope starts with a smooth sine squared function, which reaches 
its maximum at three cycles, and afterwards the temporal envelope is constant at its 
maximum value, until the last three cycles when it ramps down, see figure [9](a). 

In figure [9](b), we plot the projection at different times during the pulse: after 3 
(when Pi = 50%), 8 (Pi = 67%), 9 (Pi = 66%) and 10 cycles (Pi = 65%). As we can see, 
after three cycles, even though the electron strongly interacts with the light field, the 
excited state remains trapped in the transverse plane, due to its ponderomotive force. 
In figure ^fc) we plot the standard deviations a x = (x 2 — (x) 2 ) and a y = (y 2 — (y) 2 ). o x 
and dy oscillate until the pulse intensity decreases, and then the electron cloud starts 
expanding linearly on both the x and y axis. This can also be seen on the last three 
snapshots of figure [9](b). 

At the end of the simulation, 30% of the wavefunction is lost in the absorbing 
boundaries, due to the long interaction. Although the wavefunction cannot expand 
in the transverse plane, it expands along the propagation direction and reaches the 
absorbing grid boundary at z = ±12.5 au. Note that in the previous simulations the 
decrease of the norm at the end of the pulse was at most 2 %. 

4. Conclusions 

We have investigated in detail the transfer of angular momentum to an atom interacting 
with an ultrashort and ultraintense light pulse carrying both orbital angular momentum 
and spin. In particular, we have performed 3D numerical simulations of a hydrogen 
atom, initially in its electron ground state, coupled to an LG pulse for different 
polarization states, carrier-envelope phases, and relative positions between the atom 
and the light vortex. We have focused on the photoionization problem by choosing the 
light carrier frequency to be uj = 1 au well above the ionization threshold of 0.5 au. 

The temporal evolution of the electron wavefunction has been projected into 
spherical harmonics obtaining perfect agreement with selection rules for light-matter 
interactions including the OAM of light. Notably, we have shown that it is possible to 
exchange more than one unit of ft of angular momentum along the propagation direction 
of the pulse fulfilling AM = £ + s where £ is the topological charge of the pulse and 
s = ±1 its circular polarization contribution. For LG pulses, the 'local' carrier-envelope 
phase is position dependent due to the LG azimuthal phase. Accordingly, we have 
checked that the photoionization probability and the transfer of angular momentum 
does not change significantly for different values of the 'global' carrier-envelope phase 
for linearly polarized LG pulses. We have also investigated the electron dynamics as 
a function of the relative position between the vortex and the hydrogen atom. For 
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small relative distances, e.g., a few au, the main dynamics are quite similar although 
small changes in the electronic transition excitations have been observed. Hence, we 
expect that far from the vortex position, the atom will not absorb OAM since it will 
interact with a plane-shape- wave field. Although the light intensity that this theoretical 
study requires is very large, its value near the pulse vortex position is very small. From 
the practical point of view, the generation of light vortices with short pulses and high 
intensities requires the use of phase holograms and achromatic set ups, as it has been 
reported recently in experiment [20J. 

To obtain physical insight into the dynamics of the electron wavefunction and the 
absorption of angular momentum, we have also computed the de Broglie-Bohm quantum 
trajectories, plotting for each scenario the dynamics of the electron in the polarization 
plane. The mean value of the position given by the trajectories coincides with the 
one obtained in the orthodox formulation of quantum mechanics, but the collection 
of trajectories gives a very different view of the process. See, for instance, the case 
of the combination of OAM and circular polarization described in figure [5j where the 
trajectories clearly rotate around the light vortex, while the mean position remains fixed. 

Finally, we would like to remark that the work presented here extends beyond 
the photoionization scenario and could be used to investigate transitions forbidden by 
standard selection rules in the electric and magnetic dipole approximation. Moreover, 
the generation of OAM laser beams with high intensity is required to excite transitions 
with large angular momentum exchange, opening the door to new applications ranging 
from high-harmonic generation with OAM [2T1 [22] in attosecond science to nuclear 
spectroscopy J23J El] . 
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Appendix A. Selection rules derivation 

In this appendix we derive the selection rules corresponding to the interaction 
Hamiltonian Hi. The derivation we present here does not differ much with the one 
in j8] , but it includes the light circular polarization which then appears explicitly in the 
selection rules. 

We start by writing the transition probability amplitude as 

tyflHMi) = -^ftk l(P ' A + A ■ PM) = - j|A£<^|r • A\fa). (A.l) 
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by using that p = — im[r, T-Lo]/h, where AE is the energy difference between states \ip[) 
and |^f). For an s-circularly polarized light beam carrying £h units of OAM, see Q, and 
assuming the dipolar approximation (A > a ) and the transverse spatial approximation 
(wq 3> a ), the atom placed in the light vortex and ignoring the temporal envelope, one 
can write (in spherical coordinates) 

V2' '" 
w 0/ 



r A 




r'^ +1 sin 



W+^e^+^e-^ + cc. (A.2) 



7T P \ \£\\ 

Thus, the angular part of the transition probability amplitude between states ipi(f) = 
U\{r) Y Li { (9, 4>) and ipf(r) = Uf(r) Y L f (#, 4>) (being Yl I {9,4>) the spherical harmonics 



with the quantization axis in the light propagation direction) obtained by inserting (A.2) 



into (A.l) gives two contributions, proportional to the integrals: 
1^,3) = Jy*^(9,4>) su/l +1 9 e^+^F L f (9, 4>)dQ, 

I 2 (£, s ) = JY*^(e,(t>) S mW +1 ee- i{i+s ^Y^(e,<p)dQ. 



(A-3) 
(A.4) 



Equation (A. 3) is the term related with e in (A.2) and is associated with the 



absorption of photons [25]. Analogously, (A.4) is related with e +luJt and is associated 



to the stimulated emission of photons. It is straightforward to see that Ii(£,s) = 
7 2 (— ^, — ^), i.e., the transitions allowed by the absorption of a (£,s) photon are the 
same as the transitions allowed by the emission of a (— l, —s) photon. 



The allowed transitions between spherical harmonics given by integral (A. 3) can 



be calculated by writing sin'^ +1 #e l( ^ +s ^ in terms of spherical harmonics and evaluating 
for which transitions I\ does not vanish (see Complements B x and Cx in [26J). Thus, 
defining AL = L{ — L\ and AM = M{ — Mi, the absorption selection rules given by 



(A. 3) are 



I AL\ < \£\ + 1 < U + U , AM = £ + s, AL + £ is odd. 



Analogously, the emission selection rules calculated from (A.4) are 

I AL\ < \£\ + 1 < L { + L f , AM = -£- s , AL + £ is odd. 



(A.5) 



(A.6) 



Appendix B. de Broglie— Bohm formulation 

In this appendix give a brief overview of the main equations of de Broglie-Bohm 
formulation of quantum mechanics P [TUl HI]. By using p = — iftV, we can rewrite 
Schrodinger equation ([!]) as: 



at 2m m 2m 

a 2 

2m 



(B.l) 



Then, by casting into it the polar form of the wavefunction, ip = Re lS / h , and separating 
real and imaginary parts, in a similar manner as done in [9|, we obtain 

dS k 2 V 2 R 1 

~dt ~ 2m~R 2m 



(VS - qAf - qV, 



(B-2) 
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Figure Bl. (colour online) Initial positions of the trajectories, on the planes (a) zx 
(b) xy and (c) yz, and (d) in space. 



OR 2 

dt 



m 



R 2 (VS - qA) 



(B-3) 



Equation (B.2) is the so-called quantum Hamilton-Jacobi equation because of its 



similarity with the (classical) Hamilton-Jacobi equation but with one additional term, 
the quantum potential, accounting for the quantum features of the system. This 
similarity suggests the definition of the particle velocity as: 



—x k [t]=v k [t] = — (VS-qA)\ 
dt m 



(x k [t],t) 



(B.4) 



Thus, (B.3) becomes a continuity equation ensuring that the trajectories distribution is 



given by R 2 (r,t) at all times. 

The initial positions of trajectories {^[io]} are distributed randomly following the 
probability density function of the hydrogen atom ground state: i? 2 (r, t ) = -\e" 2r ^°, 
see figure Bl. Then, after solving ( |B.l ) (or, alternatively, (B.2) and ( |B.3[ )) we find the 
quantum trajectories, {#&[*]}, (time evolution of the positions) by integrating ( B.4[ ): 



x k [t] = / v k [t]dt. 

JtQ 



(B.5) 
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